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(57) Abstract: The invention relates to a method of reconstruct- 
ing an MR image (14) from MR signals (1 , 2, 3, 4) that are ac- 
quired in parallel by a plurality of receiving coils (5, 6, 7, 8) with 
incomplete sampling of the spatial frequency space, the recon- 
structed MR image (14) being calculated iteratively as the solu- 
tion of a system of linear equations. To provide a more highly 
developed method for parallel SENSE imaging that is an im- 
provement on the prior art in respect of the computing time re- 
quired for image reconstruction, the invention proposes that the 
iteration process (15) begin with a starting image (9) which ap- 
proximates to the MR image to be reconstructed (14), in such 
a way that undersampling artifacts are suppressed in the inter- 
mediate solutions of the system of equations that are obtained 
during the course of the itemdon process. 
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Method of reconstructing an MR image 



The invention relates to a method of reconstructing an MR image from MR 
signals that are acquired in parallel by a plxirality of receiving coils with incomplete sampling 
of the spatial frequency space, the reconstructed MR image being calcxilated iteratively as the 
solution of a system of linear equations. 
5 The invention also relates to an MR apparatus for performing the method and 

to a computer program for an apparatus of this kind. 

In MR imaging, the location of the magnetization of nuclei within the 
examination volume is usually performed by means of temporally variable, spatially 
inhomogeneous magnetic fields (magnetic field gradients). The MR signals used for 

1 0 reconstructing the images are recorded, in the time domain, in the form of a voltage that is 
induced in a receiving coil arrangement surrounding the examination volume under the 
influence of a suitable sequence of magnetic-field gradient ptilses and high-frequency pulses. 
The actual reconstruction of the images is then performed, for example, by Fourier 
transformation of the time-related signals. The sampling of the spatial frequency space (what 

15 is called the k-space), by which space the volumetric region to be imaged, or field of view 
(FOV), and the image resolution are in turn defined, is determined by the sequence that is 
preset for the magnetic-field gradient pulses and high-frequency pulses. The requirements 
that haye to be met by the FOV and the image resolution preset the number of phase 
encoding steps in the sampling of the spatial frequency space and hence the duration of the 

20 imaging sequence. A major problem with MR imaging arises as a direct result of this. It often 
takes an undesirably long time to acquire an image of the desired resolution of the complete 
FOV. This is particularly true when it is not only two-dimensional MR images that are to be 
acquired but also three (or more) dimensional ones such as, for example, ones having 
additional spectral dimensions. 

25 A large number of fiirther technical developments in the field of MR imaging 

are aimed at shortening the image acquisition time. Developments to the apparatus that make 
it possible for the magnetic field gradients to be switched as fast as possible have today 
reached the limits of what is feasible in technical terms and of what is physiologically 
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acceptable for the patient. However, for a large number of applications the acquisition tinoies 
are still too long. 

There appears to be some prospect of overcoming the technical and 
physiological limits mentioned to which conventional MR imaging is subject with the 
5 parallel methods that have become known recently, such as, for example, the so-called 

SENSE (sensitivity encoding) technique. The finding on which this latter method is based is 
that the spatial sensitivity profile of the receiving coil arrangement used impresses on the MR 
signals spatial information that can be used in the reconstruction of the images. When a 
plurality of separate receiving coils each having a different spatial sensitivity profile are used 

10 in parallel, i.e. simiiltaneously, then, by combining the MR signals that are acquired in each 
case from an incomplete sampling of the spatial frequency space, the acquisition time for an 
image can be reduced in comparison with the conventional methods by a factor that, in good 
cases, is equal to the number of receiving coils used. 

Because, in SENSE imaging, particular importance is attached to ttie spatial 

15 sensitivity profiles of the receiving coils used, the image reconstruction from the MR signals 
that are acquired in parallel by incomplete sampling of the spatial frequency space cannot 
take place, as is usual in the conventional MR methods, solely by Fourier transformation. In 
fact, special signal processing techniques are required for the reconstruction of the images in 
SENSE imaging. 

20 For cases where the sampling of the spatial frequency space is Cartesian, there 

are known reconstruction methods for SENSE imaging in which a conventional Fourier 
transformation is first carried out for each receiving coil and this produces image data sets, 
affected by undersampling artifacts, whose number corresponds to the mmiber of receiving 
coils used. In a second step, these image data sets are then combined to give a definitive MR 

25 image, with the knowledge that exists of the spatial sensitivity profiles of the receiving coils 
being used to allow the undersampling artifacts to be removed. 

If, in the SENSE imaging, the sampling of the spatial frequency space takes 
place in a non-Cartesian fashion, such as, say, radially or spirally, the reconstruction of the 
images proves to be substantially more complicated. PrUssmann has proposed a method that 

30 is particularly effective in numerical terms in which the reconstructed MR image is 

calculated as the solution of a system of linear equations (see Klaas P. Prussman et al.: 
"Advances in Sensitivity Encoding with Arbitrary K-Space Trajectories" in Magnetic 
Resonance in Medicine, Vol. 46, pages 638 to 651, 2001). Because of the size of the system 
of linear equations, it is solved, with an acceptable amount of computing work, iteratively by 
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what is termed a conjugate gradient (CG) method. A disadvantage of the method proposed in 
the article mentioned is that what the iteration process begins with as a starting image is a 
"zero image", i.e. an entirely black image, which is a poor approximation of the definitive 
MR image. What this means is that, taking this starting image as its point of departure, the 
5 known iteration process converges comparatively slowly. Correspondingly lengthy 
computing times are required to reconstruct the definitive MR image. 

Against this background, it is an object of the present invention to provide a 
method of reconstruction for SENSE imaging that is an improvement on the prior art in 
respect of the computing time required for image reconstraction. 

1 0 Starting firom a method of the kind defined in the opening paragraph, this 

object is achieved by virtue of the fact that the iteration process begins with a starting image 
which approximates to the MR image to be reconstructed, in such a way that undersampling 
artifacts are suppressed in the intermediate solutions of the system of equations that are 
obtained during the course of the iteration process, 

15 It is found that, when a "zero image" is used as a starting image for the 

iteration, the intermediate image that exists after the first step of the iteration shows 
pronounced undersampling artifacts. A comparatively large number of iterations are then 
required before such undersampling artifacts have been entirely removed and the definitive 
MR image is available as a solution. In accordance with the invention, a considerably smaller 

20 number of iterations are needed for reconstructing the MR image if what is selected as a 

starting image is an image that is a better approximation of the MR image to be reconstmcted 
than the zero image that is usually used. What in fact this ensures is that undersampling 
artifacts are suppressed in the intermediate images that are calculated during the course of the 
iteration process. The advantageous possibility exists of presetting accordingly the starting 

25 image with which the iteration process begins. What the invention achieves in this way is that 
the computing work required for the reconstruction of the definitive MR ims^e is greatly 
reduced in comparison with the prior art. 

In an advantageous variant of the method according to the invention, the 
starting image is generated by reconstructed restricted sub-regions of the starting image firom 

30 those MR signals that are received by receiving coils that are of high respective sensitivities 
in the sub-regions concerned. Hence the starting image is assembled fi:om part-images whose 
number corresponds to the nxmiber of receiving coils used, each of the part-images being 
reconstructed directly firom MR signals associated with individual receiving coils. The sub- 
regions are selected in this case in line with the geometry of the spatial sensitivity profiles of 
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the receiving coils concerned. Particularly fast convergence of the iterative reconstruction 
process, for moderate reduction factors, is obtained by this procedure. 

Cases are conceivable in which, due to the degree of undersampling, certain 
sub-regions of the starting image cannot be reconstructed in the manner described above 
5 directly from the MR signals acquired by individual receiving coils, without the starting 
image suffering to an excessively high degree from undersampling artifacts. However, even 
when there is a high degree of imdersampling (a high reduction factor), the starting image can 
be directiy reconstructed wdth a small number of undersampling artifacts at least in individual 
sub-regions. The possibility then exists, by taking the image information for tiiese sub- 

10 regions as a basis, of estimating the undersampling artifacts that can be expected in the sub- 
regions of the starting image that are still missing. Use can .be made in this case of the 
knowledge that exists of the sampling pattern that is used in acquiring the images, i.e. the 
knowledge in particular of what is called the point-spread function. The missing image 
information in the starting image can then be synthesized from the MR signals by correcting 

15 the image reconstructed directiy from the MR signals in the missing regions in line with the 
estimate of the undersampling artifacts that will occur in these regions. After this, tiie image 
as a whole will have a minimum of imdersampling artifacts. 

Alternatively, the starting image can be reconstructed at a resolution that is 
reduced in comparison with the definitive MR image. The desired resolution is then only 

20 obtained in the further course of the reconstruction process. 

In cases in which either a plurality of MR images having FOVs that are 
adjacent in space are generated, to allow a volumetric image to be generated from slice 
images, for example, or a series of MR images spread over time are generated, in the imaging 
of dynamic processes, for example, a possibility that is open is, in accordance with the 

25 invention, to use as the starting image an MR image that already exists and that is sitnilar to 
the MR image to be reconstructed at the time. The MR image that has already been 
reconstmcted does not contain any xmdersampling artifacts whatever, which means that, by 
using such an image as the starting image, particularly fast convergence is ensured for the 
method of reconstruction according to the invention. Hence, what may be used as a starting 

30 image is, for example, an MR image that has already been reconstructed and whose image 
plane is adjacent to the MR image that has yet to be reconstructed. In acquiring a series of 
MR images spread over time, what can be used as a starting image is an MR image that was 
acquired previously as part of the series and that has already been reconstructed. It is also 
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conceivable for the starting image to be generated by interpolation or extrapolation from 

existing MR images that have already been reconstructed. 

The method of reconstruction outlined above is suitable above all when the 

sampling of the spatial frequency space for the acquisition of the MR signals takes place 
5 radial or spirally or in some other non-Cartesian fashion. It is, however, also conceivable for 

the method according to the invention to be used in cases where Cartesian sampling of the 

spatial frequency space is used. 

An MR apparatus suitable for performing the method according to the 

invention has a main field coil for generating a homogeneous static magnetic field in an 
1 0 examination volume, a plurality of gradient coils for generating magnetic field gradients in . 

the examination volume, at least one emitting coil for generating high-frequency fields in the 

examination volume, a pliirality of receiving coils for the parallel acquisition of MR signals 

from the examination volume, and a central control unit for operating the gradient coils and 

the emitting coil, plus a reconstruction and display unit for processing and showing the MR' 
1 5 signals. The method described above can be performed on the MR apparatus according to the 

invention by means of a suitable programmed control means for the reconstruction and 

display unit. 

The method according to the invention may be made available to the users of 
pieces of MR apparatus in the form of a corresponding computer program. The computer 
20 program may be stored on suitable data carriers, such as CD-ROMs or floppy disks, for 

example, or it may be downloaded to the reconstruction and display unit of the MR apparatus 
over the internet. 

These and other aspects of the invention are apparent from and will be 
elucidated with reference to the embodiments described hereinafter. 

25 

In the drawings: 

Fig. 1 is a flowchart of a first variant of the method according to the invention. 
Fig. 2 shows a second variant of the method according to the invention. 
30 Fig. 3 shows an MR apparatus according to the invention. 

Figure 1 shows the procedure followed in the method according to the 
invention of reconstructing an MR image. The point of departure for the reconstruction is 
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four MR-signal data sets 1, 2, 3, 4. The MR signals have been acquired in parallel, each by 
incomplete radial sampling of the spatial frequency space. The arrangement of the receiving 
coils 5, 6, 7 and 8 relative to the circular imaging region that is to be examined (the FOV) is 
shown in the figure. A starting image 9 is generated from the MR-signal data sets 1, 2, 3 and 
5 4 by reconstructing restricted sub-regions 10, 1 1, 12, 13 of the starting image 9, which 
restricted sub-regions are acquired by the receiving coils 8, 5, 6 and 7, which are of high 
respective sensitivities in these sub-regions. To allow a definite MR image 14 to be 
reconstructed, an iteration process that is known from the prior art is begun with the starting 
image that has been generated in this way, the process being indicated in the figure by the 
10 symbol 15. The iteration process 15 produces the definitive MR image 14 as the solution of a 
system of linear equations. The fact that the starting image 9 is generated in the maimer 
described above from the MR signals ensures that undersampling artifacts are suppressed in 
the intermediate solutions of the system of equations that are obtained during the course of 
the iteration process 1 5. 

15 In the variant of the method according to the invention that is shown in figure 

2, the MR-signal data sets 1, 2, 3 and 4 are acquired with a particularly high degree of 
undersampling, i.e. with a high reduction factor. What this means is that, in a sub-region 16, 
the starting image caimot be reconstructed, in accordance with the invention, directly from 
the MR signals acquired by the individual receiving coils 5, 6, 7 and 8, However, the missing 

20 image information within the starting image 9 can be synthesized in the sub-region 16 from 
the MR signals in their entirety. This is done, as described above, by estimating from the sub- 
regions 10, 11, 12 and 13 the undersampling artifacts that can be expected in the sub-region 
16, thus enabling the starting image 9 to be corrected in such a way that it contains as few 
imdersampling artifacts as possible. 

25 Figure 3 shows an MR apparatus in the form of a block diagram. The MR 

apparatus comprises a main field coil 17 for generating a homogeneous static magnetic field 
in an examination volume in which a patient 1 8 is situated. The MR apparatus also has 
gradient coils 19, 20 and 21 for generating magnetic field gradients in different directions in 
space within the examination volume. The pattern followed by the magnetic field gradients 

30 with time and in space within the examination volxmie is controlled by means of a central 

control unit 22, which is connected to the gradient coils 19, 20 and 21 via a gradient amplifier 
23. Also forming part of the MR apparatus shown is an emitting coil 24 for generating high- 
frequency fields in the examination volume. For the parallel acquisition of MR signals from 
the examination volume, use is made of three receiving coils 25, 26 and 27, these latter being 
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so arranged that they have different spatial sensitivity profiles within the examination 
volume, thus enabling them to be used for SENSE imaging. The emitting coil 24 is connected 
to the control unit 22 via an emitting unit 28. The MR signals acquired by the receiving coils 
25, 26 and 27 are demodulated and amplified by means of a receiving unit 29 and fed to a 
5 reconstruction and display unit 30. The MR signals processed by the reconstruction and 

display unit 30 can be shown by means of a screen 3 L The reconstruction and display unit 30 
has a suitable programmed control means to allow the method described above to be 
performed. 
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1 . A method of reconstructing an MR image (14) from MR signals (1 , 2, 3, 4) 

that are acquired in parallel by a plurality of receiving coils (5, 6, 7, 8) with incomplete 
sampling of the spatial frequency space, the reconstructed MR image (14) being calculated 
iteratively as the solution of a system of linear equations, characterized in that the iteration 
5 process (15) begins with a starting image (9) which approximates to the MR image to be 
reconstructed (14), in such a way that undersampling artifacts are suppressed in the 
intermediate solutions of the system of linear equations that are obtained during the course of 
the iteration process (15). 

10 2, A method as claimed in claim 1, characterized in that the starting image (9) is 

generated by reconstructing restricted sub-regions (10, 1 1, 12, 13) of the starting image (9) 
from those MR signals (1, 2, 3, 4) that are acquired by receiving coils (5, 6, 7, 8) that have 
high respective sensitivities in these sub-regions. 

15 3. A method as claimed in claim 2, characterized in that, in those sub-regions 

(16) of the starting image (9) that cannot be reconstructed directly from MR signals (1, 2, 3, 
4) acquired by individual receiving coils (5, 6, 7, 8) owing to the degree of imdersampling, 
the missing image information is synthesized from the MR signals in their entirety. 

20 4 . A method as claimed in claim 1 , characterized in that the starting image is 

reconstructed from the MR signals at a resolution that is reduced in comparison with the 
definitive MR image. 

5. A method as claimed in claim 1 , characterized in that what is used as the 

25 starting image is an MR image that has already been reconstructed and that is similar to the 
MR image that has yet to be reconstructed. 

6. A method as claimed in claim 5, characterized in that, in the reconstruction of 
a plurality of MR slice images, the starting image is generated from at least one MR image 
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that has akeady been reconstructed and that has an image plane adjacent to the MR image 
that has yet to be reconstructed. 

7. A method as claimed in claim 5, characterized in that, in reconstructing a 

5 series of MR images spread over time, what is used as a starting image is an MR image that 
has already been reconstructed and that was acquired at a time prior to the MR image that has 
yet to be reconstructed. 

8. A method as claimed in any of the foregoing claims, characterized in that, in 
10 the acquisition of the MR signals (1, 2, 3, 4), the sampling of the spatial frequency space 

takes place radially or spirally or in any other non-Cartesian fashion. 

9. An MR apparatus having a main field coil (17) for generating a homogeneous 
static magnetic field in an examination volimie, a plurality of gradient coils (1 9, 20, 21) for 

1 5 generating magnetic field gradients in the examination volume, at least one emitting coil (24) 
for generating high-frequency fields in the examination volume, a plurality of receiving coils 
(25, 26, 27) for the parallel acquisition of MR signals fi'om the examination volume, and a 
cential control unit (22) for operating the gradient coils (19, 20, 21) and the emitting coil 
(24), plus a reconstruction and display imit (30) for processing and showing the MR signals, 

20 characterized in that the reconstruction and display unit (30) has a programmed control 
means that operates by the method claimed in any of claims 1 to 8. 



25 



10. A computer program for an MR apparatus as claimed in claim 9, characterized 

in that a method as claimed in any of claims 1 to 7 is implemented on the reconstruction and 
display unit of the MR apparatus by the computer program. 
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